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91.   Introduction. 

Meteorolo/;,ists  have  during  raany  decades  tried  to  construct 
laboratory  models  of  the  --arth' s  atmospheric  circulations. 
Since  IJorld  '.;ar  II  these  attempts  have  met  with  notable  success 
in  a  number  of  diffei-ent  places  and  much  useful  information 
besrin'i  on  v.'erther  processes  and  general  circulation  patterns 
has  been  obtained.   Of  course  the  atmosphere  cannot  be  modeled 
exactly  in  the  labor-story  in  the;  sense  of  making  all  dimension- 
less  scale  parameters  have  the  same  value  for  the  model  as  for 
the  atm.osphere.   Nevertheless,  it  has  been  shov.Ti  by  Pultz  and 
pthers  [2,  3s  k,    5,  7]  that  if  attention  is  focused  on  particu- 
lar aspects  of   the  atmospheric  motions  certain  scale  parameters 
are  very  much  more  important  than  others.   For  example,  the 
Rossby  number,  which  measures  the  ratio  of  the  inortial  force 
to  the  Coriolis  force,  is  by  far  the  most  important  scale 
parameter  for  the  modeling  of  tlie  large  scale  general  circula- 
tion patterns. 

Two  of  the  laboratory  models  which  permit  the  duplication 
of  the  Rossby  nu'.iber  make  use  of  rotating  liquids  in  cylindrical 
vessels.   Both  have  hor-n   very  successful  in  duplicating  the 
hemispheric  quasi-stationary  pattern  of  cyclones  and  anticyclones 
which  is  a  feature  of  mid-latitucle  vjeather  maps.   In  the  first 
of  these,  the  'heating'  model,  the   flow  taices  place  in  the 
annular  rejicn  betv.-een  tuo  cylindrical  surfaces  kept  at  different 
temperatures,  the  colder  surface  being  nearer  the  axis.   Because 
of  differential  volume  changes  due  i:o  the  heat  transfer  between 


the  inner  and  outer  cylinders,  a  therraal  vjind  is  created  v;hich 
produces  a  continuous  wind  shear  and  vorticity  change  in  the 
vertical  direction.   In  the  second  or  'polar  front'  iiiodel,  two 
iromlscible  liquids  of  very  sli5;htly  different  densities  are 
rotated  in  the  cylindrical  pan.   The  velocity  and  vorticity  are 
constant  along  verticals  in  each  liquid  separc'tely  but  have  a 
jump  discontinuity  at  the  interface  between  the  two  liquids. 
All  of  the  wind  shear  is  conco^.trabcd  in  the  interface.   The 
energy  needed  to  drive  thiC  fl:Lid  in  this  model  comes  from 
changes  in  the  pan  rotation  speed. 

Both  of  the  iiodels  described  above  are  baroclli:iic,  i.e,, 
the  surft-ces  of  constant  pressure  do  not  coincide  with  those  of 
constant  volume,  or  to  put  it  more  meaningfully,  the  vorticity 
is  not  constant  along  verticals.   The  models  are  therefore 
capable  of  displaying  some  of  tiie  most  complex  features  of 
atmospheric  circulations,  such  as  the  deepening  and  filling  of 
cyclonic  systems,  as  well  as  the  simplest  features  like  relative 
wave  soeed  and  amplitude.   The  'heating'  model  is  somevjhat 
closer  to  the  Larth's  atmosobcre  in  its  mode  of  operation  and 
is  also  somewhat  easier  to  manipulate  experimentally.   On  the 
other  hand,  the  'polar  front'  model  seemed  to  me  to  offer  cer- 
tain advantages  froi}  the  point  of  view  of  a  thorough  mathemati- 
cgl  treatment,   I  suspect  very  strongly,  however,  that  the 
mathematical  techniques  developed  here  will  also  be  of  use  in 
studying  the  'heating'  model. 

The  sequence  of  events  which  occurs  in  the  'polar  front' 
model  has  been  recorded  by  P'ultz  [5,  6,  9]  but  I  should  like  to 
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recapitulate  it  here  briefly  because  of  its  sin;nif icance  for 
the  subsequent  oabheraatlcal  dcvolo-oraent .   The  apparatus  is 
simplicity  it;::clf^  conrxsting  of  a  cylindrical  -_^las3  pan  Mounted 
axially  on  a  table  v/hose  rotation  rate  can  bo  controlled  over 
wide  rallies.   The  liquid  used  for  the  bottoiii  layer  is  a  monochlo- 
robf.nzeiif-toluene     solution  whose  density  can  be  adjusted  to 
a  value  just  sli/^htly  greater  than  that  of  the  water-deter/vent 
mixture  ;jhica  is  used  for  the  upper  layer.   Because  of  the  very 
small  density  difference  (about  u,5  percent  in  most  of  the 
experiments)  the  loi-;er  liquid  has  very  nearly  zero  v/eicht.   The 
shape  of  the  interface  betvieen  the  two  liquids  is  therefore  very 
sensitive  to  small  accelei-ations.   From  rest  the  liquids  are 
gradually  brou'^ht  to  a  stead^.^  rotation.  At   the  slow  speeds 
which  are  used  the  free  surface  is  scarcely  curved  but  the 
interface  is  a  rather  steeo  paraboloid  dished  upward.   The 
pressure  .gradient  force  at  the  interface  is  in  equilibriam  ^^;lth 
the  centripetal  acceleration  and  tl'io  quite  sraall  vjoight  of  the 
lower  liquid.   The  slopes  ai-e  therefore  of  the  order  of  1  rather 
than  l/lOO,  as  they  are  in  the  atmosphere  along  frontal  surfaces. 
Unless  the  lov/er  lavei'  is  quite  deep  and  the  an^ailar  speed 
small,  the  denser  liquid  will  ordinarily  be  forced  away  from 
the  axis  into  an  annular  rinj'  at  the  lower  edge  of  the  cylinder 
with  the  interface  touching  the  bottom  of  the  pan  in  a  circular 
trace . 

If  radial  equilibrium  with  essentially  solid  body  rota- 
tions in  both,  layers  is  achieved  at  a  relatively  low  speed  and 
the  pan  is  subsequently  speeded  up  by  20  or  30  percent,  the 
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motions  in  the  two  layors  take  on  quite  a  different  aspect. 
Because  the  lov;er  liquid  adjut;ts  itself  more  readily  to  the  new 
rotation,  wherees  the  u,)oer  liquid  seems  not  to  sense  the  change 
In  speed,  a  velocity  difference  is  created  bctvjeen  the  tvjo 
layers.   The  first  effect  of  this  is  to  force  the  ring  of  dense 
liquid  even  farther  awa-'-  f ro n  the  axisj  the  second  effect  of  the 
shear  is  to  produce  waves,  by  mcc'ianisras  which  are  still  not  too 
well  understood,  uoon  tie  interf "ice,   Althougi'i  the  viave  raotions 
extend  throughout  bof'  liquids,  tliO'j  are  to  an  extent  corapensat- 
ing,  since  the  free  surface  changes  relatively  little  from  its 
Initial  paraboloidal  shape  of  slight  curvature.   The  number  of 
frontal  waves  and  the  soced  vrith  which  they  move  relative  to  the 
mass  motion  of  the  two  layers  depends  critically  upon  the  shear 
that  has  been  introduced.   The  waves  theviselves  are  remarkably 
steady,  but  as  the  secondare'  frictional  circ^ilations  in  the  two 
layers  gradually  erase  the  volocicy  difference,  tiiere  are  rather 
abrupt  transitions  into  nev/  modes  of  oscillation  vjith  an  evolu- 
tion tovjard  fewer,  slower  moving  vjaves  of  longer  vjave  length. 
The  late  stages,  as  might  be  ex;~iected,  are  inclined  to  be  con- 
plicated  and  to  a  large  extent  mipredictable , 

The  type  of  flovj  described  above  is  called  abnormal 
because  it  is  essential!^'-  the  opposite  of  mr.t:ions  encoimtered 
on  the  Larth.   The  so-called  normal  flow  will  be  obtained  in  the 
model  if  the  initial  radial  equilibrium  is  acldeved  at  a  rela- 
tively high  rotatiou  r-te  and  t'le  pan  is  subsequentlv  slovjed 
dovm  by  20  or  3C  p  rccnt.   ■i?h_e  first  effect  of  the  sloiung  down 
of  the  lower  layer  -s  to  cause  it  to  collect  in  the  center 


forming  a  dorae  of  dense  liquid  on  the  axis.   The  shape  of  the 
dome  in  the  early  sta^^es  is  a  slightly  flattened  paraboloid 
dished  downwards.   The  second  effect  of  the  velocity  shear  is 
a'-ain  to  produce  lAjaves  upon  the  frontal  ouarface,  the  number  and 
V78ve  speed  a.'-ain  defending  upon  the  ai-aount  of  shear  which  has 
been  introduced.   /  feature  of  the  motions  is  the  formation  of 
a  rather  pronounced  westerl-"-  .ict  or  zone  of  high  speed  flovj  over 
the  region  of  greatest  slooe  on  tb,e  interfac"^,  which  occurs 
usually  alon-T  the  pt3ripheral  ed-je  of  the  dome.   This  same  phe- 
nomenon occurs^  though  not  so  noticcabl^',  in  tl}e  abnormal  flow 
regime.   Such  beliavior  is  not  unexpected  in  viex-j  of  the  other 
evidences  of  a  quasi-  eostrophic  adjustment  of  the  velocity  and 
pressure  fields.   The  evolution  of  the  motion  tov/ard  fewer  and 
longer  x-javes  is  the  saiae  as  beiore  but  there  is  an  even  rreater 
tendency  in  the  norrml  flow  regime  for  nonlincrr  effects  to 
come  into  play.   These  produce  a  distortion  of  the  x;ave  trou;:hs 
with  a  steepening  of  tlte  vjestern  slooe  and  Tlatteninc^,  of  the 
eastern  slopes,  producin--  the  fa:ailiar  classic  cold  front  - 
Warm  front  picture.   Eventuall"-  occlusion  of  the  xjaves  occurs 
and  the  motion  degenerates  into  turbulence. 

The  sequence  of  events  as  outlined  above  is  very  suggestive 
for  the  mathematical  development.   The  fact  that  equilibrium 
flows  exist  and  are  stcible  suggests  that  xie   should  look  for 
exact  solutions  of  the  equations  of  motion  and  boundary  condi- 
tions which  depend  only  on  the  radial  coordinate  and  the  height, 
This  is  carried  out  in  ivrtlcle  3.   i-.lthough  tlie  equations  are 
to  an  extent  under-determined  all  the  solutions  have 


the  property  that  the  interface  nay  be  dished  either  up  or 
down  depending  upon  the  relative  Magnitude  of  the  transverse 
velocities  in  the  two  layers  and  their  density  ratio.   The  most 
useful  equilibriura  solution  seems  to  be  the  one  discovered  by 
Helmholtz  in  1888  [1]  corresponding,  to  solid  body  rotations  in 
the  tv;o  layers. 

Because  so  many  features  of  the  actual  motions  are 
explained  at  least  partially  by  the  solid  body  equilibrium  solu- 
tions it  seems  reasonable  to  look  for  perturbation  solutions 
which  are  in  some  sense  close  to  the  solid  body  rotations.   The 
simplest  raetliod  of  insuring  this  property  is  to  take  the  first 
variation  of  the  full  differential  equEitions  and  boundary  condi- 
tions and  insert  in  then  for  the  unvaried  quantities  the  solid 
rotation  solutions.   The  pertarbc,tion  equations  wLich  I  obtain 
in  this  way  bear  a  rather  striKin;^  resenblance  to  the  equations 
of  motion  vjhich  are  used  in  dynamical  rieteorolo^y .   Although  the 
perturbation  equations  are  linear,  they  are  still  partial  differ- 
ential equ,':tions  ±n   three  independent  variables  r,  9,  t  and  are 
six  in  number,   A  frontal  assault  on  the  equ^tions  to  determine 
solutions  satisfvin"  ;iven  initial  cojiditlons  seems  therefore  to 
be  out  of  the  question. 

A  prominent  fe  ture  of  trie  observed  wa^-e  motloas,  however, 
was  their  rem.^rkable  stabilitv,  ir.dicetin":  the  possibilitj'-  that 
progressing  wave  type  soluoions  might  e;:ist.   I,  therefore,  look 
next  for  a  factorization  of  the  j.- rt^rbation  qu;-.ntities  into  a 
form  factor  dependent  only  on  the  radii-.s  jiultiolied  by  a  periodic 
function  of  a  linear  corabinatlon  of  an  ,le  v^.riable  and  the  time. 
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The  partial  differential  equations  are  thereby  reduced  to 
coupled  ordinary  differential  equations.   One  further  feature 
of  the  observed  wrves  has  still  to  be  incorporated  into  the 
mathematical  theory:   the  compensation  mechanism  which  permits 
waves  of  lar-e  ai-:plitude  on  the  interface  but  leaves  the  free 
surface  of  the  upper  liquid  relatively  unperturbed,  although 
the  motions  in  the  upper  liquid  are  far  from  negli-tible .   On  the 
other  hand  simply  setting  the  perturbation  of  the  free  surface 
equal  to  zero  is  equivalent  to  makin^  the  perturbation  of  the 
interface  also  equal  to  zero  urileos  the  i.ave  speeds  are  restrict- 
ed in  a  certain  way  by  sett i 13  one  parameter  in  the  theory  equal 
to  zero.   The  real  .iustif ication  for  this  whole  procedure  is  to 
be  found  in  the  fact  that  when  the  density  ratio  of  tlie  two 
liquids  is  nearly  unity  the  type  of  perturbation  which  actually 
occurs  corresponds  to  a  singular  perturbation  of  the  differential 
equations  in  the  small  parameter  e  =  1-P,  vihere  p  is  the  ratio 
of  the  density  in  the  upper  liquid  to  that  in  the  lovjer. 

After  these  preliminaries  what  is  left  is  mathemst  ics. 
The  result  of  tl'se  assu'rptions  above  and  som.e  eli-,;inations  is  an 
ordinary  difre-.-'cntia]_  equation  of  second  order  for  the  radial 
dependence  of  the  height  of  the  interface.   After  some  not  quite 
obvious  changes  of  v;.rieble  t:.e  equation  can  be  brought  to  the 
standard  form,  of  the  h-rpergeometrio  differential  equation.   It 
turns  out  that  a  solution  satisfying  the  boundai'y  conditions  at 
the  outer  wall  of  tije  cylinder  and  regular  on  the  axis  can  exist 
for  a  ;:iven  number  of  waves  around  the  periphery  of  the  cylinder 
only  if  the  angular  speeds  in  the  two  layers  are  in  a  certain 


ratio  to  each  other.   Tiie  actual  numerical  value  of  the  ratio 
is  found  from  a  rather  coirrollcatec'  transcendental  equation. 
The  results  are  suriiinarized  in  Tables  I-III.   The  explanation  of 
their  calculation  is  rjiven  in  Part  II  of  this  report  along  uith 
graphs  of  the  solutions  and  a  detailed  analysis  of  their 
correspondence  with  the  observed  results. 

^ ^ •   i-atheraatical  formulation  of  the  lorobleiri . 

In  order  that  the  results  can  be  iMrnediately  apolied  to 
an  analysis  of  the  J'ultz  experiments,  I  shall  derive  the  equa- 
tion of  motion  in  terms  of  cylindrical  oolar  coordiiiates  r,  9,  z 
with  corresponding'  velocit'^'-  components  u,  v,  v/.   The  velocities 
will  in  all  cases  be  referred  to  the  nonrotatin^  laboratory 
reference  frarie.   In  a  subsequent  report  it  xjill  be  shovm  how 
a  very  similar  development  can  be  carried  out  for  the  exterior 
of  a  herriispherc  vjith  rravity  directed  radiall"  ini-jard  and  rota- 
tions occurring  about  tlie  polsir  axis.   The  ti/o  situations  are  in 
fact  not  greatly  dissimilar  since  \ie    can  ima.'r,ine  a  continuous 
mapping  or  deformation  of  th.e  physical  space  ivhich  XNfould  leave 
the  axis  of  rotation  ir.variant  but  vjould  carry  the  plane  hori- 
zontal sui'faces  of  tiie  cylindrical  flox-j  over  into  the  concentric 
spherical  horizontal  surfaces  of  the  liemisplier ical  flow.   This 
basic  physical  sinilaritv  accounts  in  larf-:e  part  for  the  rele- 
vance of  the  Fultz  experiments  to  the  study  of  atmospheric 
circulations.   Fluid  dynamical  IJ.hcnesses  are  also  significant, 
of  course. 
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The  pressure  will  be  denoted  by  p  and  the  constant  density 
by  p  (its  reciprocal,  the  specific  volurae  a,  ulll  also  be  used 
occasiono.lly ) .   In  order  to  dXotin:_;uish  the  two  layers  of 
liquid  I  shall  use  a  prime  on  all  quantities  referring  to  the 
upper,  less  dense  liquid  while  tbe  absence  of  a  prime  refers  to 
the  lox^jer,  denser  liquid. 

The  boundaries  of  the  domain  are  determined  by  the  bottom 
surface  z  =  0,  the  cylinder  surface  r  =  R,  the  free  surface  of 
the  upper  liquid  z  =  ^(r,S,t),  and  the  interface  z  =  4'(^j^jt) 
separating  the  two  liquids.   On  these  boundin;;:  surfaces  we  have 
the  folloivinj  dynamic  and  hineraatic  conditions  v;hich  must  be 
satisfied  at  all  times; 

(2.1)  p«  =  0   and   w'  =  (L  ■:-  u'^  +  r~"^v''(lJ,,    on   z  =  (f 

(2.2)  w  =  0    on    z  --^  0 

(2.3)  u  =  (J   and   u '  =  0    on   r   -  II 

{2.1].)  P  =  P'   and    w  ~  ucj' ,  -  r  '  v4c.  =  W  -  n'4,  -  r   v'cJ)q 

on   z  =  4 

Finally  we  shall  im.pono,  for  reasons  of  mathematical  simplicity, 
a  kind  of  syimnetry  on  the  motions  to  be  considered  by  requiring 
that 

(2.5)  u  =  0   and    u'  ^  0    on   r  =  C   . 

The  equations  describiri^r  the  incompressibility  of  the  txjo 
liquids  are 
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(2.6)  u   +  ^~   u  +  r'  v^  +  w  =  0 

(2.7)  u>  +  r"  u'  +  r"  V'  ^-  v."  =  0 

r  9    z 

The  conservation  of  morientun  is  sir.iil&rly  described  by  the 
equations 

(2.8)  U-,  +  uu.   +  r~  vu„  +  wu   -  r"  v   +  ap  =0 

t     r        0     z  ^T 

(2.9)  V,  +  uv^,  i-  r"  vv^  +  wv  +  p~  uv  +  ar~  p^  =  0 

t     r        9     z  9 

(2.10)  p  +  CD   --  0 

^    -  z 

(2.11)  u'  -;-  u'u'  +  r   v'u'  -;-  vj'u'  -  r~  v'''"  +  a'^'  =  0 

t      r         9      z  '■1"' 

(2.12)  V'  +  utv'  +  r~  v'v'  +  vj'v'  +  r'^^u'v'  +  a'p'^^p'  .=  0 

(2.13)  g  +  a'pi  =  0   . 

z 

(The  hydrostatic  approximation  that  vertical  accelerations  are 
snail  compared  to  the  ^Travity  acceleration  has  been  used  in 
equations  (2.10)  and  (2.13).) 

Although  the  differential  equations  and  boundary  conditions 
above  co-niiDletely  describe  the  physical  problem  subject  to  appro- 
priate initial  conditions  they  are  because  of  their  number  and 
complexity  unsuitable  for  the  fiirther"  discussion  of  the  "oroblem. 
Our  object  will  be  to  elirainate  so:':e  of  the  dependent  variables 
Qnd  incorporate  the  boundsry  conditions  on  the  bottom,  the  inter- 
face and  the  free  surface  explicitly  Into  a  new  set  of  equations. 
This  object  can  be  achieved  by  two  different  methods.   The  first, 
described  in  more  detail  in  Appendix  A,  employs  the  additional 
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assumption,  amply  justified  by  the  observations,  that  all  the 
shear  of  the  horizontal  velocity  components  is  concentrated  in 
the  interface.   This  method  has  the  advantage  that  it  leads  to  a 
set  of  six  nonlinear  differential  equations  ^^7ith  independent 
variables  r,  9  and  t.   These  equations  would  be  useful  for  num- 
erical computation  of  the  solution  for  given  initial  conditions 
and  for  analyzing  sbability  of  various  ty:>es  of  motion. 

The  second  method,  to  be  described  in  Article  q.,  has  bhe 
advantages  for  our  present  purposes  that  it  does  not  require  any 
additional  assumptions  and  proceeds  directly  to  a  set  of  simpli- 
fied linear  equations, 

§ 3 .   Equilibrium  sol ution 3  of  solid  body  ro tati o n  type . 

It  is  readily  verified  that  equations  (2,1)- (2.13)  allow 
exact  solutions  in  which  the  radial  and  vertical  components  of 
velocity  vanish  identically  and  in  x-«Thich  the  remaining  dependent 
variables  are  independent  of  the  tii:ie  and  angle  variables.   For 
these  solutions  the  pressiire  gradient  in  each  liquid  is  in  exact 
balance  with  the  centripetal  acceleration  directed  raoially  and 
the  gravity  acceleration  directed  vertically.   The  solutions  are 
underdetermlned  to  the  extent  that  the  transverse  velocity  dis- 
tribution can  be  prescribed  as  an  arbitr.'iry  function  of  the 
radius.   Thus  there  exist,  for  example,  double  layer  flox-;s  of 
the  rtankine  vortex  t^^^pe  in  which  the  velocity  increases  linearly 
from  zero  at  the  center  to  a  maximum  at  some  prescribed  distance 
and  then  decreases  like  the  inverse  first  povrer  of  the  radius 
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out  to  the  edre  ox"  the  cYlinc-er.   All  of  the  conceivably  signi- 
ficant flows,  hovjever,  shai'c  their  must  important  proosrties 
with  tlie  siiiiplest  type  of  solution,  that  corresponding  to  solid 
body  rotations  of  different  speeds  in  the  tuo  layers.   The  rest 
of  the  discussion  will  be  liraited  to  such  linear  velocity 
profiles. 

It  is  readily  verified  that  corresponding  to  the  prescribed 
transverse  velocity  distributions 

(3.1)         V  =  ior   ,    V'  =  ,  i'r   , 

where  'a^-   and  ui'  are  constant  angular  speeds,  we  have  the  height 
and  pressure  formulas 

1,2  2 


(3.2)  e|(r)  -  gq;^  +  i  of    r 

(3.3)  g4(r)  =  g(L  +  ~   2fiL-r„P 


:l.,,l 


where  ^  =  i!:(^')  ^'^^^  *i   ~  4('  )  ^^e  prescribed  constants,  and 

O.k)  P'  (r*,z)  =  rjP'  (|-z) 

(3.5)  p(r,z)  ■--  gp'(3--cl>)  +  Sp(4-2)   . 

Some  interesting  properties  of  the  double  layer  solutions 
are  easily  deduced.   Preferring  to  ecm;'.tion  (3.2)  \-ie   note  that 
the  height  of  the  free  surface  is  independent  of  the  densities 
in  either  layer,  that  it  is  a  paraboloid  dished  upward  and  that, 
if  account  is  taken  of  the  exnerimental  magnitude  of  the  angular 
speed  (iJ  ,    the  curvature  is  very  small.   Prom  equation  (3.3)  we 
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note  that  the  interface  is  likeuiss  a  paraboloid,  but  in 
contrast  to  the  free  surface,  its  slope  and  curvature  depend 
strongly  on  the  densities  and  on   the  an-^.ular  speeds  in  the  two 
layers.   In  the  experiiiients  the  density  difference  is  small  and 
therefore  if  the  an;:7alar  speeds  are  not  also  nearly  the  sarae  the 
interface  vjill  be  sci-ongly  curved,   nor  cover,  since  the  angular 
speed  in  the  bottom  layer  cari  bo  either  r;reater  or  less  than 
that  in  th^  top  layer,  the  paraboloid  can  be  dished  either  up 
or  down.   These  properties  of  the  solution  are  in  complete 
accord  with  the  exoerlraental  observations. 

For  convenience  later  we  Introduce  the  thicl-rness  4'  -   i!l~4 
of  the  u;iper  layer  ,'-iven  by  the  forraula 

(3.6)  ^4,(r)  .  ,4:  -^^^i^^t^^^r^      . 


§l\..      Perturbations  fron  solid  rotation. 

In  order  to  find  a  simpler  set  of  linear  equations  which 
describe  motions  close  in  sorae  sense  to  the  solid  body  rotations 
of  the  double  layer  we  may  take  the  first  variation  of  the  com- 
plete equations  and  bonndsry  conditions  of  §2  and  afterwards 
substitute  for  the  unvaried  quantities,  wiierevcr  they  occur,  the 
solid  body  solutions  of  the  last  section,   ^.'e  denote  a  perturb- 
ation qua;itity  by  ore  fixing  a  "5''  to  the  symbol  for  tl.e  corres- 
ponding unvaried  dependent  variable.   It  is  convenient  also  to 
introduce  the  notation  for  the  linear  operators 
D  =  o'/3t  +  ^u  ?/'^Q   and  D'  =  ?/^t   s  '.jJ    ^/de.   The  result  of  the 
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operations  described  above  on  the  equations  (2.1)- (2. 13)  leads 
to  the  linear  equations 

(,(;.. 1)  6p»  =  0   and  5vj'  =  D'o^  +  ^^ou'     on   z  =  1^  +  Scf 

([|..2)  Sw  =  0    on  2  =  0 

(lj..3)  5u  =  0   and  5u»  -  0    on   r  =  R 

(I).,!;)  p  +  6n  =  p»  +  op'  arid 

ov/  -   4  ,5u   -   i'/64ci   ~   '->■'■'    -   4   Su'    -    i.O'oc|^        on      z   =  4   +  o4 
r  y  r  t? 

(lj..5)  5u  -  0        and        5u»    =0  on       r   =  0 

([i..6)        6u      +  r'   ou  +  r'   5v,,   +   ow     =0 
r  9  z 

ik'l)        6u'    +  r'^^ou'    +  r'-^ov'    +   5w'    =   0 
r  y  z 

([|.,8)  D5u  -    2u.i5v   +  a5p  =   0 

{l\.'9)  D5v   +  2uj5u   +  ar"-^5p      =  0 

(k-lO)  a5p  =   0 

ZJ 

(l|..ll)         D'5u'    -    2t'"ov'    +   a'5p'  =0 

([;.12)         D'ov'    +   2u''5u'    4    a'r"-'-o-o^,   =   0 
(1;.13)  a'6p^         =  0 

Some   eleraentary  operations   lead  at   once   to  the   results 
ik'lk)  a'5p'    =   r6|  =  g(64    +54') 
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(i;.l5)       a6p  =  gpfef  -  54)  +  8^4  =  g(o4  +  P54') 

where  p  =  p'/p  -  a/a' .   Also  the  velocity  perturbations  5u,  6v, 
5u'  and  ov'  are  found  to  be  independent  of  z  so  that  the  vertical 
velocity  perturbations  can  be  found  explicitly  in  terras  of  the 
horizontal  velocities: 

('|.16)     ow  -  -z(ou  +  r~  5u  +  r"  ovr) 

(ij..l7)    OW  ^  D'Ofl"  -:-  ((r-z)(Ga'  +  r"-'-ou'  +  r"-'"uv')  +  cf^uu' 

Finally  the  equations  satisfied  by  the  six  perturbation  quanti- 
ties ou,  ov,  54j  5u' J  uv '  aiid  o4'  ?re 

('|.l8)  D6u  -  2u.6v  +  r:(o4  +  Po4')^  =  0 

([j..l9)  D5v  +  2w6u  +  .-2;r"^(o4  +  p64')^  =  0 

([j_.20)  D64  +  4  ou  +  4(Qu  +  r'-'-ou  +  r"-^6v^}  =  0 

(,'4.21)  D'ou'  -  2u."ov'  +  £-(54  +  5f|'')p  =  ^^ 

([{..22)  D'uv'  +  2..:i'6u'  +  n-r"--(64  +  ^4')@  =  0 

(L.23)  D'54'  +  4-'ou'  +  4' (6u.'  +  r~-'"5u'  +  r'-'-Sv')  =  0 

Surmac.rizinf;  the  resul.ts  just  obtained,  we   may  say  that  the 
assuraptions  of  ner;ll2;iblo  vertical  acceleration  (hydrostatic 
balance)  and  snail  deviations  from  rotationally  symiiietric  motion 
lead  at  once  to  tne  co''iplete  vertical  hcraogeneity  in  each  layer 
and  the  concentre tlon  of  all  velocity  shear  in  the  interface. 
Vie  note  also  the  characteristic  property  that  in  the  lower  layer 
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ascending  raotion  is  associated  x-:ith  convergence  of  the  velocity 
field  and  descending  motion  witli  divergence;  the  magnitude  of 
the  vertical  velocity  is  greatest  at  the  interface.   finally  it 
is  worth  pointing  out  the  striking  similarity  betxjeen  the 
equations  ([}..lO )  -  ('1..23 )  and  the  equations  of  dynamical 
meteorology. 

§ S .   Progressing  wa ve  so 1 utions . 

The  experiments  of  ?ultz  indiccte  ijlainly  tht.t  steady  i^jave 
motions  of  the  ti;o  layers  are  possible  or  rather  would  be  so  if 
friction  could  be  entirely  eliminatec'.   I  therefore  look  for 
such  progressing  waves  travellijig  with  a  uniform  angular  speed 
k:.   The  solutions,  if  they  exist,  will  then  depend  periodically 
on  the  variable  ^  =  Q  -  Kt.   If  there  ere  an  integral  number  m 
of  wave  crests  around  tlie  circu/iforence  we  may  write 

(5.1)  -[ou,Cu'5  (r,S,t)  =  {u,U'i(r)  sin  my 

(5.2)  {ov,5v' ,b4,o(j)'i  (r,e,t)  =  {VjV' ,H,  H' Wr  )  cos  me   . 

From  now  on  we  shall  indicate  derivatives  with  respect  to 
r  by  a  dot  over  the  quantity.  VJith  the  substitutions  above  the 
final  equations  of  the  last  section  become 

(5.3)  r^-(W-  -  u^)^  +  2..'V  =  g(h  +  pii') 
(5.[|.)       2'ujU  +  m(K  -  u!)V  =  giftt-^'^di  +  pH'  ) 

{S»s>)  ra(«.  -  u')ll  +  4'J  +  4(U  -•.-  r"^U  -  mr'-'-V)  =  0 
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(5.6)  ra(!i    -    w''  )U'    +   2i.'V'    =   g(H  +   H«  ) 

(5.7)  2.<Ui    +  r,i(iL    -   .X'')v'    =   ,;,riir"-^(H  +   H' ) 

(5.8)  iri(ic    -    i'J  )IV    +  c|)'U'    +  4' ('J'    +  r"   U'    -  wr" V)    =  0      . 
If  we  I'lake   the   abbreviations 

(5.9)  lv  =  m^dt    -    ij')^   -  kic^ 

(5.10)  ^'    =  m^(^    -    u')^   -   W^ 
then,    providing  O"  ^  0 

(5.11)  U  =  Q-r\[{rc   -      ')(H  +   pi')    -    2c:r~^(n  +   PH')] 

(5.12)  V   =   --'r^[2'o(H  +   pn«  )    -   r(i^(K:   -    io)r'-^{'ii  +  pH' )  ] 

(5.13)  U  =   -T"^ia[  (Vv.  -  uo)  (I-I  +  pli'  j    -   2u.r"-'-  (H  +  pfl'  )    +  2o^r"^(PI  +  pH'  )  ] 

These  expressions  may  be  substituted  into  equation  {S'5)    to 
obtain  the  differential  equation 

(5.1^)   4[(H+PK')  +  r"^(H  +  PN')  -  m^r'^dl  +  pE-  )] 

4-  4[(n+Pl*I')  -  2"(K..  -  ..■)~-^r~^(H  +  PH'  )]  +  s'\e  =   0  . 

It  is  an  observed  fact  in  the  Fultz  experliraents  that  when 
the  density  difference  e  =  1-p  is  snail,  the  perturbations  of 
the  free  surface  are  negligible  compared  to  the  oerturbations  of 
the  interface.   This  prorerty  is  not  evident  fron  the  equations 
but  it  is  physically  plausible  since  the  Interface  would  be  very 
sensitive  to  s/jall  accelerations  when  the  lower  liouid  had 
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nearly  zero  weight,   Refcrr3.n<'3  to  equations  i$.G)    and  (5.7)  we 
see  that  when  H  +  H'  is  set  equal  to  zero  the  velocities  U'  and 
V  satisfy  a  pair  of  homogeneous  algebraic  equations.   For  non- 
trivial  solutions  to  exist  tlie  deterrp.inant  of  coefficients,  i.e. 
0-'  defined  in  equation  (5.10),  must  vanish.   V.'e  thus  have  the 
relations 

(5.15)  ra(k  -  u;'  )  =  2sa;< 

(5.16)  u»  =  -sv 

vjhere  s  =  -i-l  or  -1, 

The  dependence  of  the  velocity  u'  on  r  is  deternined  after 
H  has  been  found  by  an  integration  of  equation  (5.8)  vjhich 
reduces  to  the  form 

(5.17)  d(r     qi'J'j  =  c:suj'iir     dr   . 

Finally  equations  (5.11 )- (5.1'!.)  all  simplify  somevjhat  when 
use  is  made  of  the  fact  that  Ii+pK'  =  (l-p)H  =  eH. 

§  6 .   The  hyp  erg  eometric  e  q  u  atio^n  . 

In  order  to  simplify  equ.ation  (5.1^)  we  make  the  substi- 
tution E(r)  =  r^'^Z(r),  so  that  Z  satisfies 


(6.1)  4Z  +  r"^[(2n+l)ci,  +  r4]Z  -!•  [r/gs  +  (m  -  j-%i;)i'""^4  ]Z  =  0   . 
Next  vjrite  the  equilitariot'a  equation  of  the  interface  in  the  form 

(6.2)  4  =  4^(1  +  kr^)   ,    k  =  (iv^  _  pco.2)/oe,.4^ 
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and  define  the  neii   independent  vai-'iable 
(6.3)  C  =--  1  -  4/4^  =  -kr-   . 

Then,  ijith  Z(r)  =  ^(?),  it  is  easily  vei-i.fied  that  ^  satisfies 
the  hyperfjeoraetric  equation 

(6.[j.)   ^(1-?)H  +  [c  -  (a+b+Dc^j^  -  ah^   --  0 

where 


(6.5) 


c  =  iTi+l  =  a+b   ,     ab  =  a  +  ?y. 


i,  =ct/2(.,2  .  p,,,„2)   ^     ^  ^  1(,.,  _  „.^) 


§ 7 .   Solution  of  the  boundary  value  problem . 

Since  c  =  m+1  is  an  integer,  it  is  knovm  from  the  general 
theory  of  the  hypergeometric  differential  equation  [8]  that  the 
approoriate  pair  of  fundanental  solutions  of  equation  (6.5)  is 
^^(5)  =  P(a,b,-o;^)  and  ^2^-)  ^  i^(a,b;l  ;1~5 ) .   The  original 
interval  for  r  is  from  r  =  0  to  r  =  R.   The  corresponding  inter- 
val  in  t,    is  therefore  from  ^^  =  0  to  £;  =  ^.-  =  -kR^  =  1  -  4p/4q« 
Of  the  two  fujidamental  solutions  only  Kiil)    'J-^  re';Tular  at  4  =  0. 

—  Til 

The  second  solution  behaves  lilce  ^   near  4  =  *^«   Since  the 
height  and  velocity  perturbations  are  required  to  be  bounded  on 
the  axis  of  the  cylinder  we  discard  the  second  solution. 

Two  cases  raust  now  be  dist int^uished  depending  on  whether 
the  interface  is  curved  upvjard  (k  >  0)  or  curved  do^^/n^•.'ard  (k<0). 
The  former  is  called  abnormal  flow  and  the  latter  normal  flovj. 
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Since  ^      >   G,  s;jn  k  =  c-n  (jj)     ~    p'v;'  ).   In  abnormal  flow 
4   <<  r))-»  v/hereas  j.n  norraal  floT/  the  o'oposlte  is  the  case, 
A     «  4  ,   The  intervals  of  inte^.Tation  in  the  tvio  cases  are: 

a)  Abnormal  flow;    -|.^^^,  |  <  ^  <  0   ,     I^p  |  -->  1 

b)  rlor.'iial  flow;     0  <  t,  <  t,^   <   1 

The  hyp er;;"eorie trie  function  ^-,(^)  is  regular  in  both  these 
intervals. 

The  boundar'''  condition  at  the  outer  wall  of  the  cylinder 
requires  that  the  radial  velocity  must  vanish.   After  the 
appropriate  transfor.iiati  ons  of  equation  (5.11)  tlie  condition 
becomes 


(7.1)  4?.  -  -XK,,  at 


4)      'op 
11 


or  by   a   well  knoxjn  oroi'^rty  of   tb.e   hyper^eoiietric   function,    cf. 
Reference    [3] 

(7.2)  abi^^,F(a+l,b+l,'    m+2;    4-.)    =   -  (rri+l )  AF(a,b;    m+1 ;    g^)       . 

Lquation  (7.2)  is  an  eigenvalue  equation  to  determine,  for 
a  given  value  of  the  integer  m,  the  ratio  of  the  speed  in  the 
lower  liquid  uo  the  S'^eed  in  tlie  upper  liquid.   The  value  of  the 
dimensionless  .Frou.de  nuniber  defined  by  the  relation 

(7.3)  yr  =   ((.;'^R^/2e34j5)^^2 

and  the  nujiber  s  =  +1  aopear  as  parameters  in  equation  (7.2) 
whose  value"  must  be  chosen  beforehand. 
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The  calculation  of  the  eigc-nvalues  and  subsequently  the 
heljhts  and  velocities  is  a  difficult  oroblera  in  itself.   The 
■procedures  which  x-jere  used  and  a  detailed  analysis  of  the  results 
will  be  ,';;iven  in  a  separate  report.   Here  I  shall  only  discuss 
the  various  cases  and  sumraarize  the  eigenvalues. 

Consider  first  the  case  of  abnormal  flow  in  which  the 
denser  liquid  is  oiled  up  near  the  cylinder  rira,   V'e  observed 
earlier  that  for  this  case  u;/y^'  >  Jp    ex  1.   From  equation  (5.l5) 
the  vjave  pro;•)a^ation  speed  is  k;/uj*  =  1  +  ~.   Thus  if  s  =  +1, 
K,'/u>'  >  1,  i.e.  the  i-javes  are  prc'^ressin^^;  the  relative  v;ave  speed 
is  highest  for  small  integer  values  of  ra  and  decreases  monotonely 
with  increasing  number  of  waves.   If  s  =  -1,  H;/^ '  <  1,  i.e.  the 
waves  are  retrograde;  for  low  integer  values  of  m  the  speed  is 
likewise  low  and  increases  with  m.   In  s.bnormal  flow  both  pro- 
gressing and  retrograde  vjaves  are  mathematically  possible.   The 
existence  of  retrograde  vjaves  in  the  laboratory  has  not  yet 
been  established  so  far  as  I  know.   Tables  I  and  II  below  give 
the  relative  wave  speeds  and  the  lower  liquid  relative  speeu  for 
various  x-jave  numbers  and  Froude  nujnbers  appropriate  to  abnormal 
flow  both  proi-ressin^;  and  retro "rade. 


ra 

u;/ujt 

h,A ' 

Fr   =   10 

Fr  •■=  l[|..l!j. 

Fr  =  17.32 

3 

1.775 

1.79I|. 

1.302 

1.667 

k 

1 .  61^0 

1,630 

1.637 

X  ,  >  U  0 

^ 

l.SkS 

1 .  330 

1.550 

l,i|.00 

6 

l.i+20 

l.i|.35 

l.li38 

1.333 

7 

1.350 

1.361 

1.363 

1,233 

« 

1.3  00 

1.308 

1.310 

1.250 

Table  I,   Progressing  waves,  abnormal  flow. 
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m 

(../.,,'.' 

yt/w' 

Fr   =   10 

Fr  =  V\..\\ 

Fr   =  17.32 

3 

2.993 

2.671 

2.516 

0.333 

k 

2. 2/^6 

2.089 

2.01[| 

0.500 

5 

2.198 

2.093 

2.0LI.3 

0.6C'0 

O 

2.21x3 

2.163 

2.127 

0.667 

7 

2.302 

2.238 

2.210 

0.717 

8 

2.359 

2.306 

2.282 

0.750 

Table  il.   lietrograde  waves,  abnormal  flow. 

In  the  case  of  norraal  flov;  the  denser  liquid  is  piled  up 
near  the  a:cis  of  the  cylinder  and  the  condition  is  that 
a)//.i)'  <  v/p  ',--■   1.   As  before  the  vjave  speed  condition  is  that 
lt/<.'->'  =  1  +  — ^.   Heuristicallj'-  it  would  seem  rather  improbable 
that  wave  speeds  greater  than  the  liquid  speed  in  both  layers 
could  exist  and  indeed  no  roots  of  the  transcendental  equation 
(7.2)  could  be  found  corresponding  to  s  =  +1,   The  normal 
progressing  viaves  correspond  to  s  =  -1.   Table  III  presents  the 
values  of  relative  wave  speed  and  relative  lower  liquid  soeed. 


m 

(O/tO' 

k/'^'' 

?1^   ^  1.1 

Pr  =1.2 

Fr  =1.3 

3 

XXX 

XXX 

XXX 

0.333 

h 

.1+33 

XXX 

x:<x 

o.5t'0 

5 

.516 

.562 

-r  -^r  -KF 

0.600 

6 

.589 

.589 

.6/^11- 

0.667 

'7 

.6  [ill 

.6!|1 

.61j.8 

0.717 

Table  III.   Progressing  waves,  normal  floTJ. 
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§8.   Appendix  A. 

If  in  the  equations  of  §2  vm   make  the  a  priori  assumption 
that  the  vertical  shear  is  concentrated  entirely  in  the  inter- 
face, i.e.  the  horizontal  velocity  coinponents  are  constant 
separately  in  each  layer  but  have  a  jump  discontinuity  at  the 
front,  then  the  conservation  of  raoraentum  equations  simplify. 
They  may  be  written 

(8.1)  Du  -  r'-'-v^  +  ap^  =  0 

(8.2)  Dv  +  r'-'-uv  +  ar'-'-p^  =  0 

(3.3)  D'u'  -  r'-'-v'"'  +  atpi  -  0 

(8.1^.)  D'v'  +  r'-'-u'V  +  a'r"  P^  =  0 

v/here  wo  have  introduced  the  nonlinear  derivative  operators 
D  =  (d/.?t  +  u  o/er  +  r""v  V--'^)  and  "'  =  (a/3t  +  u'  3/ar 
+  r'-'-v'  "/OG). 

The  equations  of  mass  conservation  in  the  tvjo  layers  can 
be  written  ir,miediately  in  the  integrated  form 

(8.5)  w  =  -z^\ 

(8.6)  w'  -  -z/y    +   4 (A'  -A)  ■"    f'^^'I'r  +  ^"■^[v]4q 

where  /\   and  /\'  stand  for  the  horizontal  velocity  diver£;ences 
in  the  two  layers  and  [f]  -  f'-f. 

At  the  interface  and  the  free  surface  the  vertical  veloci- 
ties can  be  v/ritten,  by  using  the  kineraatic  boundary  conditions, 
as 
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(8.7)  l^;  (z  =  4)  -  d4 

(8.8)  w'  (z  =  f )  =  D'lF  =  D'4  +  i3'4' 

Ttien  it  follovjs  froui  these  equations  and  equations  (8.5)  and 
(8.6)  applied  at  z  =  (j>  and  z  =  Ip  respectively  that 

(8.9)  d4  +  4(u  +  r"-^u  +  r"-''v,^)  -  0 

(8.10)  D'4'  +  4'(u'  H-  r~-^u'  -!-  i'""''v')  =  0 

Integration  of  the  hydrostatic  pressure  equations  gives  as 
before 

(8.11)  a'p'  ^-  g(4  +  4'  -  z) 

(8.12)  o,p  =  3(4  -i-  p4'  -    z) 

which  raay  be    substituted   in  equations    (8 .1 )- (G.li.)    to    complete 
the    systeri. 

It    is    soraetimes   useful   to    consider   the   vorticities    in  the 
tv/o   layers,    defined   by 

(8.13)  11    "^  '^r   "*"  ^"''^   "   ^''"   '■'■© 
(8.11^)  Sil=  v^   -1-  r"^v'    -   r'^-ut^ 

Then   it    is    easily   deduced   fron  the   preceding  results   that   the 
potential  vorticities  il /4   ^nd    i ,:.'/4'    ^r^    covjscrved   in  the 
raotion,    i.e. 

(8.15)  Diri/4)    =  0 

(8.16)  "'  (:yv'4')  =  0 
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